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Abstract 

We analyse the elastic energy of an intercalated smectic where 
orientationally ordered polymers with an average orientation varying 
from layer to layer are intercalated between smectic planes. The lowest 
order terms in the coupling between polymer director and smectic 
layer curvature are added to the smectic elastic energy. Integration 
over the smectic degrees of freedom leaves an effective polymer twist 
energy that has to be included into the total polymer elastic energy 
leading to a fluctuational renormalization of the intercalated polymer 
twist modulus. If the polymers are chiral this in its turn leads to a 
renormalization of the cholesteric pitch. 







Recent elucidation of the structure of DNA - cationic lipid complexes 
JD, || has brought forth quite a few unsuspected features of this macromolec- 
ular aggregate. It appears that cationic lipids in the complex retain their 
preferred packing characterised by a multilamellar 3-D smectic order while 
the oppositely charged DNA gets intercalated in between the lipid bilayer 
smectic planes. The intercalated DNA appears to be packed with a 2-D 
smectic order of relatively small domain sizes that are probably coupled to 
the order in the neighboring intercalated DNA layers || . The careful! X-ray 
diffraction studies leave no ambiguity as to the fact that both components 
of the complex - cationic lipids as well as DNA are ordered. 

The ordering tendencies giving rise to this complex aggregate structure 
are due partly to the fairly well understood interactions between lipid bi- 
layers in aqueous solutions where van der Waals attraction competes with 
electrostatic and hydration forces, augmented by entropic repulsion forces 
originating in elastic fluctuations of the lipid bilayers constrained to a mul- 
tilamellar stack ||. Similarly interactions between DNA molecules in the 
bulk that have only recently come under closer experimental as well as the- 
oretical scrutiny, appear to be dominated by repulsive forces of electrostatic 
as well as hydration origin, this time too augmented through entropic mech- 
anism very similar to the one operating in multilamellar lipid systems [II] . 
The interaction between the two constituents of the complex are probably 
dominated by the electrostatic attraction between DNA and cationic lipids 
possibly modified by elastic shape fluctuations of the DNAs intercalated be- 
tween positively charged layers of lipids and possibly by the forces mediated 
by the lipid bilayer elasticity due to local deformations induced by the close 
proximity of intercalated DNA. More work is certainly needed to asses the 
relative importance of all these mechanisms in bringing about the stability 
of the DNA - cationic lipid complex. 

Compared to the phases existing in the bulk the lipid subphase does not 
appear to be substantially modified. It has the same structural geometry 
as the one found with other lipids in the bulk. DNA is in this respect very 
much modified. At effective interhelical spacings found in the DNA - CL 
complex 0, DNA in the bulk would be either in the line hexatic phase 
or within the cholesteric phase ||. Very little of this bulk order persists 
in intercalated DNA that is forced into effectively 2-D layers intercalated 
between lipid bilayers. The positional order does not change qualitatively 
if we consider only DNA intercalated within a single layer. They are both 
short ranged ||. The orientational order is changed more drastically if the 
state of affairs in the bulk and in the DNA-CL complex are compared. It is 
nevertheless the apparent total absence of the cholesteric order in the complex 
that seems to us the most baffling. Apart from very tentative statements || 
that cholesteric structures of extremely large pitch ( ~ mm) can sometimes 
be detected in the complex the chiral nature of the DNA molecule makes no 
imprint on the structure of this macromolecular aggregate. 

It is our goal in this contribution to investigate the interaction between 
orientational order of intercalated polymers and the smectic degrees of free- 
dom of intercalating lipid bilayers. We propose a simple theory of the effect 
that the coupling between polymer (i.e. DNA) orientational ordering within 
the intercalated layers and smectic order of these layers can have on the 
effective twist elastic constant of the polymer layers. If in addition the in- 
tercalated polymers are chiral, this theory for the first time introduces a 
comprehensive mechanism for coupling between smectic and cholesteric de- 
grees of freedom leading to a fluctuational renormalization of the cholesteric 
pitch of the polymer subsystem. 

Model We will consider a simplified model of an intercalated smectic 



the director of the polymers n(n; p) within an intercalated layer is a constant. 
We introduced n as the height index of the layer (layers are assumed to have 
the average positions at z n = n x d where d is the average layer - layer 
separation) while p — (x, y) is the transverse radius vector. We have thus 
effectively limited ourselves to a mean-field approximation within a single 
layer. 

We will construct an elastic free energy of this system, assuming a general 
dependence of the polymer director on the position i.e. height index n, of the 
layer. Since the polymer orientational order can interact with the curvature 
energy of each layer we first have to construct all the scalar invariants that 
can be composed with n and the 2nd fundamental form of a single smectic 
layer. 

If one defines the 2nd fundamental form of Gauss of the n-th layer with lo- 
cal displacement described within the Monge parametrization (x, y, ( n (x, y)) 
as 

d 2 ( n (p) 

where the indices i, k can have values i, k = 1,2 with x\ — x and X2 = y, 
then the lowest order scalar invariants which can be built from the director 
n and the tensor Eq.[TJ of that layer are three || 
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Only two of these invariants are linearly independent if we ignore the terms 
containing Gaussian curvature. The curvature elastic energy of n-th surface, 
assuming that within each layer the director of the polymers is a constant 
n(n; p) = n(n), can be written as 

H n (K ik ;n) = \K C J ' (TrK lk (n; p)) 2 d 2 p + 

+ \a J K u (n;p)K ik (n;p)n l (n)n k (n)d 2 p + 

+ \b J K ik (n; p)Ki m (n; p)ni(n)n k (n)ni(n)n m (n)d 2 p 

(3) 

In a stack of layers the total elastic energy is composed of curvature elastic 
energy Eq.||] and the deformation energy due to smectic dilations - compres- 
sions of the layers in the transverse direction, characterised by a smectic 
compressibility modulus B 

H = j H n (K lk ; n)dn + \B JJ (^^j dn d 2 p (4) 

The formal limits of this expression are well known and will not be discussed 
here. Introducing Fourier transform of the local displacement in the direc- 
tions p with wave vector Q we obtain for the total smectic elastic energy the 
expression 
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where we have defined the operator H(n, ri; Q) and used the shorthand J2q = 

I d 2 Q with S being the area of the layer. The elastic energy expression 

Eq. [| represents the final formalisation of our model. If the intercalated 
polymers are chiral we have to add to Eq. |5] the standard cholesteric elastic 
energy which depends in the lowest order only on the derivatives of the 
polymer director with respect to the stack index ||. 

Free energy of smectic fluctuations We now proceed by integrating out the 
smectic fluctuations from the free energy defined through elastic Hamiltonian 
Eq.|| and thus obtaining an effective intercalated polymer free energy that 
will depend only on the director field of the polymers. We start by setting 

0(n;Q) =aQ 2 (Qn(n)) 2 + 6(Qn(n)) 4 (6) 
and writing out explicitely the operator we introduced above Eq. [5] 

H(n, ri; Q) = + $ K <& + Q)l 6 ( n ~ ( 7 ) 

which allows us to express the free energy corresponding to smectic elastic 
fluctuations as 

JF(n(n)) = - kT In 5 = - kT In (Yl Q j ' ■ . j V( n (Q) exp{-(3 #)) = 
kT 

= — ^hiDefflfan'iQ). (8) 
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In order to get the part of the free energy that depends explicitely on the 
polymer director field, it is convenient to resort to the following representa- 
tion of the free energy Eq.|8] 

kT 
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where T$ is the part of the free energy that does not contain n(n) explicitely. 
The Green function Q^n, ri; Q) entering the above equation can be obtained 

as pn 

{^PB^p + P K cQ 4 + vP<t>(n; Q)) G„(n, ri; Q) = S(n - ri). (10) 
Introducing now Qo(n, ri; Q) = Go(\n — n'\; Q) = Q^ =Q (n, ri; Q) in the form 

G (n - ri; Q) = |\/JV 2 exp S^Q 2 \n - ri\ (11) 

we can expand Q^(n,n'; Q) perturbatively up to the first order in <p(n; Q) 
thus obtaining the free energy to the second order in this quantity 



Since the zero order Green function is obviously of short range we can expand 
the expression for the free energy Eq.[L2| for a slowly varying <p(n; Q) field in a 
standard fashion, obtaining the following approximate form of the free energy 



^(n(n)) = Jo + iE k(0; Q) / dn 0(n; Q) - 

Q L J 

r dtQo(t\ Q)) J dn (j) 2 (n; Q) + 



+ f ( fdtt 2 g 2 (t;Q)) [dn 



dn 



(13) 



The summation over the Fourier space intends also integration over the dif- 
ferent directions of the Q vector. Denoting this orientational integration with 
<. . .>(jj it is easy to see that <4>(n; Q)>cl> as well as <4> 2 {n; Q)>cl> do not 
depend on the orientational angles at all and are thus independent of the 
director field n(n). The dependence on the director field remains only in the 



derivative terms, i.e. terms of the form ( 9 ^q^ ) • These terms contain 



d(Qn(n)) dn(n 



dn 



Q- 



dn 



Qri(n) = Q (n(n) x fi(n)) 
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where H(n) is the vector of the "angular velocity" of rotation of the director 
from layer to layer. If the average normal to the layers is in the z direction 
than O(n) = (0,0,Q z (n)), also h(n) is within each layer and in direction 
prependicular to n(n) with magnitude |ii(n)| = Q z (n). As can be easily 
seen, the only terms that survive the integration <. . .>qj and still depend 
on the director field are those depending quadratically on |n(n)|. From the 
free energy Eq.^ these terms can be obtained in the form 

^•(n(n)) = *b(a, 6)+|^ J^Q'dQ (J dt t 2 g 2 (t; Q)) ((a + bf 

' ' '(15) 

where ^(a, b) is the part of the free energy that after the <. . integration 
does not depend explicitely on the director any more. Evaluating the last 
integral over Q and taking into account the fact that the minimal value of Q 
is set by the domain size, assumed to be a square of side /, while the maximal 
value is set by the molecular dimension a, we obtain finally 
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Clearly for this geometry K st is the additional twist elastic constant of in- 
tercalated polymers, stemming from the smectic interactions between layers 
of orientationally ordered polymer molecules. Thus we see that in an inter- 
calated smectic the smectic interactions tend to renormalize the twist elastic 
modulus of intercalated oriented polymers to 
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This renormalization is of purely fluctuational origin. If the polymers are in 
addition chiral, the renormalization of the twist modulus obviously leads 
to unwinding of the cholesteric pitch of the intercalated polymers to a new 
equilibrium value of 
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Discussion The mean - field model introduced above led in a very straight- 
forward way to a renormalization of the twist modulus of intercalated poly- 
mers. The effect itself is a very intuitive one. If there exists a deformational 
"easy axis" within each smectic layer, which according to Eq.[5] represents a 
deformational wave whose direction is perpendicular to n(n) within a single 
layer, the smectic compressibility term would tend to twist the directors of 
the neighboring layers towards a colinear position. This would introduce a 
coupling term (n(n) — n(n + l)) 2 in the free energy of two neighboring layers. 
The continuum version of this effect would lead exactly to Eq.|T6|. 

The coupling constants a and b, c.f. Eq|| between the orientational or- 
dering of intercalated polymers and effective elastic properties of the layers 
depend in general on the orientational order parameter of the polymer chain 
S, the elastic modulus of the chains defined as kTC p where C p is the per- 
sistence length and the polymer surface density p. The orientational order 
parameter of the polymer chains is defined through the 2-D orientational 
tensor <Tjfc of the intercalated chains as a a = |(1 ±5), where a is the index of 



the two eigenvalues |TTJ. For S close to 1, i.e. close to complete orientational 
order where all the chains point in the same direction, the scaling form for a 
and b should be 

a,b~kTS^, (19) 

where £± is the separation between the chains perpendicular to their long 
axis. This is the form appropriate for our assumption of complete ordering 
of chains within each smectic layer. 

We can now asses the magnitude of the contribution of smectic modes to 
the twist elastic modulus of the smectic layers. Assuming the scaling Eq. [19] 
for constants a and b we obtain 

Kst ~ {kT)~ l (-^-) kTS^) (20) 
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times an unknown numerical constant persumably on the order of 10. In the 
units used in Eq. [| -j^- is dimensionless. Since the dominant interactions 
determining B are electrostatic attractions between DNA and cationic lipid 
headgroups the smectic modulus has to be quite large while -f- is on the 



order of 10, the modification of the twist modulus implied by Eq. |T7| thus 
has to be enormous. If the intercalated polymers are chiral it wouTa thus 
come as no surprise if the effective cholesteric pitch surviving in this system 
would be orders of magnitude larger than in the pure polymer system |J. 

We have not included the possible direct modification of the chiral in- 
teractions by smectic fluctuations in this analysis, an effect that would act 
in the direction opposite to the smectic fluctuation renormalization of the 
polymer twist modulus. As this effect would tend to make the effective pitch 
smaller, which apparently has never been observed in this system, we assume 



pitch in such a system, if the polymers themselves are chiral. We propose 
this as the primary reason why no cholesteric structures have been observed 
in the DNA-CL system. 
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Figures Fig.l. A schematic representation of a part of a polymer inter- 
calated smectic system. The smectic layers are not shown explicitely. The 
direction of the intercalated polymers changes from layer to layer. It has 
been arbitrarily assumed to change by j between two neighboring layers. 
The smectic interactions tend to orient the neighboring polymer layers in a 
parallel direction. 
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